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We study the properties of generalized striped phases of doped cuprate planar quantum anti- 
ferromagnets. We invoke an effective, spatially anisotropic, non-linear sigma model in two space 
dimensions. Our theoretical predictions are in quantitative agreement with recent experiments in 
La2-iSr a: Cu04 with < x < 0.0f8. We focus on (i) the magnetic correlation length, (ii) the 
staggered magnetization at T = and (iii) the Neel temperature, as functions of doping, using 
parameters determined previously and independently for this system. These results support the 
proposal that the low doping (antiferromagnetic) phase of the cuprates has a striped configuration. 
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There is no direct evidence for periodically structured 
striped phases in doped antiferromagnetically ordered 
cuprates, such as La 2 _ x Sr x Cu04 (with x < 0.02). In- 
deed, recent neutron scattering measurements [jjj show 
that the incommensurate magnetic peaks which are char- 
acteristic of dynamic stripe formation at higher hole con- 
centrations x, seem to disappear below x of order 0.05, 
well above the antiferromagnetic regime. Yet the growing 
theoretical literature Q supporting the tendency toward 
microscopic phase separation of holes in these strongly 
correlated quasi-two-dimensional materials, as well as the 
experimental observation of stripe phenomena in 

many related systems, raises the strong possibility that 
the doped antiferromagnetic cuprates will also be char- 
acterized by such stripes. Moreover, there is indirect evi- 
dence of linear, or striped, features in the magnetic struc- 
ture of these materials, including the finite size scaling 
properties j^] of the Neel temperature and uniform mag- 
netic susceptibility, and the successful interpretation |||7]] 
of muon spin resonance (/iSR) and nuclear quadrupole 
resonance (NQR) experiments within models that pre- 
sume a striped structure. 

Those latter theories are nominally based on a static 
periodic array of stripes, which separate antiferromag- 
netic slabs, or ladders, at most weakly coupled to one 
another across the stripes. Since such regular arrays are 
not observed in the neutron scattering experiments, we 
want to consider a broader class of striped structures. 
These may include arrays with varying separation be- 
tween neighboring stripes, as suggested Q by the neu- 
tron scattering lineshapes in the related doped nickelates, 
or dynamic behavior of the stripes, including amplitude 
fluctuations, rigid translations, or the meanderings pro- 
posed Q by Zaanen and coworkers. We also note that a 
magnetic phase domain, rather than anti-phase bound- 
ary p0[ ], which locally only suppresses, rather than re- 
verses the antiferromagnetic order parameter, would give 
at best only a weak incommensurate scattering, even for 
a periodic array. It is then our purpose here to develop a 
reliable effective field theory which will predict the exper- 
imentally observable behavior while making a minimum 



of assumptions about the details of the underlying struc- 
ture, beyond the demand that it corresponds on average 
to periodic one-dimensional weakening of the antiferro- 
magnetic exchange between spins. 

Encouraged and guided by the remarkable success 
of Chakravarty, Halperin and Nelson (jl]J in describ- 
ing undoped La2Cu04 by just such an approach, we in- 
troduce a suitably anisotropic non-linear sigma model 
parametrized so as to reproduce the long wavelength be- 
havior of the doped material in the antiferromagnetically 
ordered phase. Whatever the origin of the average pe- 
riodic modulation of exchange, the corresponding long 
wavelength magnetic excitations from the antiferromag- 
netic ground state, which determine the thermodynam- 
ics, will have the same character as those of a static ar- 
ray of stripes with the same period. They are gapless 
(Goldstone modes, from spin rotational symmetry), with 
a dispersion relation for the frequency squared which is 
analytic in wave vector, giving a linear dispersion with 
anisotropic "spin wave" velocity, where the principal axes 
of the velocity tensor must be parallel and perpendicular 
to the stripes (which lie along a crystal axis j| ) . Such be- 
havior is described in the continuum limit by the action 
of an effective anisotropic non-linear sigma model, 
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where h is a unit vector field. The symbols have been 
chosen to suggest the continuum limit of an underlying 
effective integer spin ff2|| Heisenberg hamiltonian on a 
square lattice with lattice constant a, and with nearest 
neighbor exchanges J x and J y . In terms of these param- 
eters the principal spin wave velocities are 
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The fundamental underlying anisotropy parameter is 



1 



then equivalently the ratio of the two exchange constants 
or of the two velocities, 
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Its value characterizes the theory, but just how it de- 
creases with increasing doping concentration x will have 
to be set later. 

We note that in restricting the subsequent analysis to 
these longest wavelength excitations we ignore both the 
non-linear dispersion of these acoustic modes and all op- 
tical branches associated with the (average) superlatticc 
in the x-direction. As an estimate of the lowest lying 
optical mode energies we can take the location of the 
lowest gap introduced by a stripe superlattice in the hole- 
free antiferromagnetic host, which occurs approximately 
at Jtt^/2/No for a superlattice of period N lattice con- 
stants in the host with exchange constant J. Coulomb 
energy costs would seem to prevent values of A*o much 
greater than 20 to 30, giving a lowest optical mode en- 
ergy greater than about 250 - 350 K. These energies lie 
above the temperatures of interest, and the neglect seems 
reasonable. 

From this point on we make use of well established 
techniques to analyze the behavior of the field theory de- 
scribed by the action (|l|) so as to make predictions for 
the physical properties of the system of interest. First 
it is useful to rewrite (|l]) more symmetrically by a di- 
mensionless rescaling of the variables: x' = (a) _1 / 4 :zA, 
y' = (a) 1 / 4 yA (A ~ 1/a is a momentum cut-off), and 
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tion (nl) becomes 



S, 



i-hkfico 




/ dr' 


[ dx' J 


Jo J 





(4) 



(aA)/S 



(5) 



(2so) Jo 

where /i takes the values x',y',T', 

g (a) = hc A/p a s =[2(l + a)/^ 1/2 

is the bare coupling constant, cq 
Jy)y / JxJ y ] 1 ^ 2 (aS)/h the spin wave velocity and p° s — 
yj J x J y S 2 the classical spin stiffness of the rescaled 
model. The original anisotropy is now hidden in the lim- 
its. We started with a problem with a finite bandwidth, 
a lower bound on length which requires us to impose a 
cutoff in the original continuum formulation. The change 
of variables introduces an anisotropy in the cutoffs. 

The (T-model action, and the spin correlations it im- 
plies, can be studied in the large N limit (TV is the number 
of components of n), where a saddle point approximation 
becomes exact 13,14], In the antiferromagnetically dis- 
ordered phase the staggered spin-spin static correlation 
function is given [jlBI by 
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where a is the component of n in the ordering direction, 



is the scaled length, rj 



2 /\/a + \fay 2 , and m is 



the inverse correlation length, given formally by imposing 
the condition that the magnitude of the field h is unity at 
each point. At zero temperature this condition becomes 
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We find that m <§C a 1 / 4 over the full range of parameters 
of interest, in which case (0) gives 
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where 



g c (a) = \/87rVa/(l + a) {in (yfa + Vl + a) 

+ Valn[(l + Vl + a)/^]}' 1 , (9) 

is the critical coupling constant of the theory: the ground 
state is disordered |16[ if m(a) > 0, or go(l) < 9c (a)- It 
decreases monotonically with increasing anisotropy 1/a, 
as the hole doping concentration x grows. Thus the sys- 
tem remains disordered at T = for a < a c , where 
a c , the critical anisotropy, is defined by the condition 
gd&c) = <?o(l)- We will see below, as this suggests, that 
as a approaches a c from above with increasing doping, 
this is also the value where the spin stiffness — and the 
three dimensional Neel temperature — vanishes. Numer- 
ically, a c « 0.047. 

For sufficiently weak anisotropy (1 > a > a c ) that 
the system orders magnetically at zero temperature, 
we can take advantage of the power of renormalization 
group techniques for systems with diverging correlation 
lengths, to understand the corrections to the classical 
limit above. As always, we introduce an explicit length 
scale, a — > Aa = e~ l , where a is the original lattice pa- 
rameter. Then the normalization condition on the field 
h at T — becomes 
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where a is the average value of the order parameter (stag- 
gered magnetization), and the second term on the right 
hand side represents the spin flip fluctuations away from 
perfect Neel order. The RG procedure can be carried 
out in various ways, including (see, e.g., [p"l| ) explicit 
integration over large momentum values and suitable re- 
scaling to restore Eq. jl0| ) to its original form. Instead 
we introduce |L8| the explicit re-scaling function Z(A): 
a = Z(A)M, and go{a) = Z(A)gR,(a), which renders 
( |l0| ) rcnormalizable, with the renormalization equation, 
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which can be integrated from the bare value go (a) at 
t = to the fully renormalized value at £ — > oo. The sta- 
ble fixed point is at gu = 0, justifying keeping only the 
lowest two powers of gn in the equation (|ll| ) , a one- loop 
approximation. The result then gives for the renormal- 
ized spin stiffness, 
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reduced from its classical value p®(a) for fixed anisotropy 
a in such a way that it vanishes, as foreseen above, at 
the critical value a — a c . 

Now we turn to prediction of and comparison with 
experimental quantities, including correlation length 
£(a,T), zero point magnetization M s (a), and Neel tem- 
perature Tjv (a.) . We know from studies of the undoped 
system that while the one-loop calculation accurately 
predicts the leading exponential dependence on inverse 
temperature of £(T), its results are not very good for the 
prefactor of that exponential or the algebraic tempera- 
ture dependent corrections, and the same will surely be 
true for the doped system, with a < 1. Therefore we 
use an interpolation formula between the exact result of 
Hasenfratz and Niedermayer JlSj] for the nonlinear sigma 
model in the neighborhood of magnetic order, and the 
result for the quantum critical regime, where £ oc T : 
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This gives excellent agreement pT[ with experiment in 
the pure case, x = 0. The experimental data Q for 
£ _1 (a;,T) have been interpreted phenomenologically, not 
according to (fL3|) , but as the sum of the pure system 
function plus a constant depending only on doping con- 
centration a;: f- x (a;,T) = f-^O.T) + f -1 (a;,0). But 
only the pure (x = 0) experimental results are reliably 
in the antiferromagnetic region we treat here. The other 
curves, with nominal values of x = 0.02, 0.03 and 0.04, 
correspond most likely to the cluster spin glass phase. 
Of course, to make a direct comparison with the experi- 
ments, yet to be done, in the region x < 0.02, we need the 
connection between the anisotropy parameter a and the 
hole concentration x. We can do this indirectly, by com- 
paring experimental properties, for example, as a func- 
tion of Neel temperature, as we will do below for the zero 
point staggered magnetization. One additional parame- 
ter is still required, though, for comparison with (13). 
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The spin stiffness ((TJ) depends through p° s on the prod- 
uct J x J y , as well as the ratio a. If we take that product, 
for example, to be independent of x, or a, then the curves 
for £ _1 (T) for doped and undoped systems will be found 
to cross, as seems to be suggested by the experiments 
@ for x = and 0.02. 

In general, the staggered magnetization M s depends 
on the short, as well as the long wavelength physics of 
the problem, so we can't use the results at T = of 



the long distance nonlinear sigma model theory directly. 
Instead, as in fll|| , we use the asymptotic long distance 
behavior of the equal time spin-spin correlation function: 
(a(x,y)a(Q,0)) -> (M s /M ) 2 as x,y ^ oo, where M 
represents perfect Neel spin alignment. To do this we 
evaluate (g) at the point where the scaled length r s is 
equal to the Josephson correlation length £j = hc/p s , 
which separates long from short wavelength scales. On 
the one hand, this is large enough for the correlation 
function to have reached its desired asymptotic behavior, 
but it is still short enough that the (exponential) decay 
due to the long wavelength fluctuations at T = + where 
(pi) holds has not yet become effective. Thus we find, 
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As usual [^TJ for quasi- two-dimensional systems, we 
make use of the long range order parameter correlations 
that have developed in the plane above the ordering tem- 
perature Tjv, to call on the validity of mean field theory 
for the weak exchange J± in the third dimension in es- 
tablishing 3-d magnetic order at a finite temperature. A 
simple physical interpretation of the result is that Tv is 
the temperature at which the energy of thermal fluctu- 
ations, fc^Tv, becomes sufficient to flip the spins in a 
region of linear dimension of the order of the 2-d corre- 
lation length, £(Tv). Since the number of spins in this 
region is proportional to (£/a) 2 , and the relative stag- 
gered magnetization in the region is given by M s /Mq, we 
estimate 



k B T N (o) w Jj 



U(T N ,q) M a 



(15) 



This expression has been used previously |11| to estimate 
J±/ks ~ 0.01 K from the experimental Tv of the pure 
material. Since M s /Mq < 1, this gives (£/a) > 10 for 
Tv > 1 K, suggesting that this mean field theory is rea- 
sonable for Tv greater than a few kelvin. 

Though the predicted staggered magnetization jl^ ) 
and the Neel temperature ( |l5| ) both depend on the 
anisotropy parameter a, whose explicit dependence on 
doping x is not established by the model, we can elimi- 
nate this dependence between the two relations, plotting 
M s versus Tv, each normalized to the corresponding un- 
doped (x — 0) value, as is done in Fig. 1. The experimen- 
tal values are those of Ref. We note, in particular, 
that the zero point magnetizations are extrapolated val- 
ues. As is explained in |J, there is a change of behavior 
around T = 30K for all values of doping. The authors 
have interpreted this as some sort of freezing of the holes, 
which may be binding to the donor impurities, estab- 
lishing static charge density waves within the stripes, or 
some other change in behavior. In any case, the reported 
values of M s are those extrapolated from the observed 
magnetization curves at temperatures above this "freez- 
ing". We emphasize that the comparison in Fig. 1 has 
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no adjustable theoretical parameters, and the agreement 
with experiment is excellent. This, and the further com- 
parisons made in 0, then supports a generalized stripe 
picture of the antiferromagnetic doped cuprates. 
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FIG. 1. Neel temperature versus staggered magnetization 
normalized to the undoped values. Line: theory; diamonds: 
experiment. 
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